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Abstract 

We review different scenarios for the motion and merging of Dirac points in two dimensional crystals. These different 
types of merging can be classified according to a winding number (a topological Berry phase) attached to each 
Dirac point. For each scenario, we calculate the Landau level spectrum and show that it can be quantitatively 
described by a semiclassical quantization rule for the constant energy areas. This quantization depends on how 
many Dirac points are enclosed by these areas. We also emphasize that different scenarios are characterized by 
different numbers of topologically protected zero energy Landau levels. 
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1. Introduction 

The spectacular properties of graphene [1] have led 
to a tremendous interest for the study of physical 
properties related to an electronic spectrum exhibiting 
Dirac points with a massless spectrum. If, in graphene, 
the two Dirac points are situated at particular points 
(K, —K) of the Brillouin zone (BZ), it is important to 
realize that in general a pair (or several pairs) of Dirac 
points may exist at any location (D, — D) in the BZ, 
and may even merge and disappear for appropriate 
change of band parameters. Moreover the exceptional 
physical properties of graphene-like systems are not 
only due to the linear dispersion relation but also to 
the two-component structure of the wave functions 
corresponding to the two atoms per unit cell. More 
precisely, each Dirac point may be characterized by a 
topological number (a Berry phase) which may take 
two values ±7r [1] . 

In this note, we review different low-energy models 
that exhibit several Dirac points and discuss different 
scenarios for the merging of these Dirac points, accord- 
ing to the topological numbers attached to each point 
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Fig. 1. Different scenarios for the merging of Dirac cones, cor- 
responding to different topological characteristics. 

(Fig. 1): a) merging of two Dirac points with opposite 
Berry phases, b) merging of two Dirac points with the 
same Berry phase, c) a more complex spectrum which 
may exhibit different types of merging, in particular 
a situation that may be called "triple merging" where 
three Dirac points merge into a single one. 

We consider quite a general structure of a two- 
component Hamiltonian for a 2D crystal with two 
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atoms, molecules or orbitals per unit cell. If inversion 
symmetry is preserved, it has the general form 



Hk 




(1) 



apart from a term /£ 1 that does not affect the spinorial 
form of the wave functions and that we omit from now 
on. The eigenenergies are ±|/fc| and the corresponding 
wave functions are given by 



V>: 



±1 



(2) 



where </>*. — arg(/&). The relative phase cj>h exhibits 
a particular topological structure: around each Dirac 
point, the circulation of <f>k along a closed path is quan- 
tized such that the topological Berry phase [2] defined 
as (J)b = \ § Vfc^fc • dk takes multiple values of it. In 
graphene, the pair of Dirac points has two opposite 
Berry phases ±7r . In bilayer graphene, the two Dirac 
points have a quadratic dispersion relation and the 
topological Berry phase is ±2tt [2,5] . Under appropriate 
conditions that we discuss later, each point ±2-7r may be 
separated into two points with a linear dispersion and 
the same Berry phase (it, it) or (— 7r, — 7r), or even four 
points with phases (it, it, it, —tv) or (— 7r, — 7t, — 7t, it) [5]. 
For the three scenarios mentioned above, we con- 
sider the spectrum without and with a magnetic field 
and show how it is constrained by topological consid- 
erations. In particular, while the maximum number of 
zero-energy Landau levels is given by wt = J^ \wi | , the 
number of such levels which are topologically protected 
is given by w p = | ^2 i w% | , which plays the role of a min- 
imal number of zero-energy Landau levels. The quan- 
tity Wi = 7^: § c Vfc(/>fc • dk is the topological winding 
number associated with each Dirac point. Moreover we 
find that the Landau level spectrum {e n } can be accu- 
rately described by the semiclassical quantization rule 



. , N n eB , 1 \wc\ 
Ac(e n ) = 2*-(n + --tf 



(3) 



where Ac(£n) is the area delimited by a closed contour 
C in reciprocal space, associated with energy e n . wc 
represents the total winding of the relative phase along 
this closed contour. 



2. (77, — 7r) — > merging 

We first consider the merging of two massless Dirac 
points with opposite Berry phases. The simplest model, 
which exhibits such a transition, is a slight modification 
of the graphene model, in which one of the transfer 
integrals is increased [6-8]: 



f k = t' + te %kai + te xka2 , (4) 

where a\ — |(e x +y/3e y ) and a^ — § (e x —y/3e y ) are 
basis vectors of the triangular Bravais lattice, in terms 
of the lattice constant a. Figure 2 shows the evolution 
of the spectrum when varying the parameter t! . When 
t' increases, the two Dirac cones approach each other 
and become anisotropic, until they merge at the time- 
reversal invariant M point when t! — 2£. At this critical 
point, the spectrum has a very peculiar structure that 
we discuss below. 
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Fig. 2. Merging of Dirac points in the toy-model of deformed 
graphene 

Starting from this toy model we have found more 
generally under which conditions, a pair of Dirac points 
may exist and merge for a general dispersion fk = 



y 



where the t mn are hopping parame- 



ters. When time-reversal symmetry is respected, Dirac 
points come always by pair (D, — D) and they neces- 
sarily merge at a time-reversal invariant point G/2 of 
the BZ, where G is a reciprocal lattice vector. In two 
dimensions, there are four different possibilities G/2 = 
(pal + qa,2)/2, where p = =bl, q = ±1 and a* are ele- 
mentary reciprocal lattice vectors. Around these points 
the Hamiltonian has the form (h — 1): 



lin 







~ ZC( iy + tt^ + A * 
2m* 





(5) 

where the parameters ra*, c, A* can be related to the 
original band parameters [9]. For example the driv- 
ing parameter of the transition A* is given by A* = 
5^m n ^mn( — l) pm+gn , where (p, q) are the coordinates 
of the merging point. Although Hamiltonian (5) has 
been derived in the vicinity of time-reversal invariant 
points, it is universal for any merging of type (tt, — 7r). 
When varying the parameter A*, this Hamiltonian de- 
scribes a topological transition between a semi- metallic 
phase with two Dirac points and an insulating phase 
with a gapped spectrum (Fig. 3). At the transition 
(A* = 0), the spectrum has a "hybrid" behavior, also 



called "semi-Dirac" : it is linear in one direction and 
quadratic in the other [7,9,10]: 



_Qx_ 
2m* 
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(6) 



We have studied in detail the Landau level spec- 
trum in the presence of a magnetic field (Fig. 4). When 
A* < 0, the two Dirac cones are well separated and 
the Landau spectrum scales as e n oc ±y/nB with a two- 
fold degeneracy corresponding to the two cones with 
winding numbers wc = =bl. When A* increases, the 
two valleys start to communicate and the two-fold de- 
generacy is progressively lifted. At the transition, the 
spectrum scales as e n oc ±[(n + 1/2) B] 2 ' 3 , as can be 
straightforwardly found from the semiclassical quanti- 
zation rule (3), where the winding number \wc\ — 
since the two Dirac points have merged and their oppo- 
site Berry phases are annihilated. Moreover the spec- 
trum is clearly separated in two different regions by 
the saddle point connecting the two Dirac cones. Be- 
low the saddle point a closed contour at a given energy 
encloses only one Dirac point and wc = ±1. Above 
the saddle point energy, such a contour encloses the 
two Dirac points and wc = 0. As shown in Fig. 4, the 
semiclassical quantization rule (3) is excellent, except 
for the vicinity of the saddle point because of the dis- 
continuous change in the area Ac and the associated 
change in the total winding number wc in Eq. (3). 
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Fig. 3. Merging of Dirac points described by Hamiltonian (5) 



In graphene, such a merging of Dirac points would 
require a modification of hopping integrals impossible 
to reach experimentally [11]. However, Dirac cones ex- 
ist in other physical systems, such as organic 2D con- 
ductors [12], optical lattices of cold atoms [13], where 
the merging transition may be more easily achieved. 
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Fig. 4. Dimensionless energy levels e' n 
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oc A*/B ' . The blue curves are the 
results from semiclassical quantization rule (3) which is excel- 
lent except near the saddle point energy which corresponds to 
the dashed blue line [9]. Here we have only shown the positive 
energy levels. The two-fold degeneracy of the zero energy Lan- 
dau level is lifted when approaching the merging transition. 

3. (77, 7r) — >> 27r merging 

We now consider the case of bilayer graphene in 
which the two layers have a translational or rotational 
mismatch (twist) with respect to the perfectly A-B 
stacked case [14]. In a first approximation (which will 
be discussed in the next paragraph), the low energy 
Hamiltonian around the K point can be written in the 
form [15] 



liq 







n 2 

2m* 

A 



(?) 



2m* 

where II = q x + iq y . When A = 0, the spectrum is 
parabolic and the associated Berry phase is 2n (— 2ty 
around the — K point). For non-zero values of A, which 
we have chosen to be real, the quadratic Dirac cone 
splits into two linear Dirac cones, though with the same 
Berry phase 7r, in contrast to the (jr, —tt) universality 
class where the relative sign in the Berry phases was 
dictated by time-reversal symmetry. Two fundamental 
differences with the previous case are noticeable. First 
the zero field spectrum stays gap less for all A, [16] 
second the spectrum in a magnetic field exhibits w p = 
2 modes of zero energy ( x4 to account for valley and 
spin degeneracy) which, in contrast to the previous 
merging, are topologically stable. 

The Landau level spectrum have been calculated in 
[15]. In Fig. 5, it is shown as a function of the dimen- 
sionless parameter A/huo c with uj c — eB/m* . At large 
A, the spectrum scales as ±vn5, as expected from 
the semiclassical quantization rule with wc — 1. When 
A is small, or well above the saddle point energy, one 
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Fig. 5. Landau levels in the model (7) of bilayer graphene 
((-7T, 7r) — >- 27r merging). The dashed lines indicate the energy 
of the saddle point, for more details see [15]. This spectrum 
has to be compared with the spectrum corresponding to the 
(-7T, — 7r) — > merging (Fig. 4). Here, the zero energy levels stay 
topologically stable. 

expects to recover the spectrum of the A-B stacked bi- 
layer, ±?ijdOc\/ri(n — 1), which in the limit of large n 
scales as ±huj c (n — 1/2) as expected from Eq. (3) with 
wc = 2. Notice that, in contrast to monolayer graphene 
where the saddle points occur at a rather high energy 
(~ 3 eV), the energy of the saddle points in twisted bi- 
layer graphene depends on the rotation angle between 
the layers and may be found in the 10 — 100 meV regime 
[20]. 



4. Multimerging 

It turns out that the spectrum in bilayer graphene 
has an interesting fine structure at low energies. Due 
to interlayer hopping between atoms which do not face 
each other, the low-energy Hamiltonian takes indeed 
the form 



/Lq 



2ra* 



2ra> 



+ vU- A 



+ vii* 



o 



(8) 



where v is a small perturbation (neglected in the previ- 
ous section) responsible for the trigonal warping of the 
isoenergy contours at low energies [5] . The parameter 
A, which we choose to be real here, mimics a strain of 
the bilayer or a displacement between the layers [17], 
with respect to perfect AB stacking (A = 0). Fig. 6 
shows the energy spectrum when A = 0. It consists of 
a central massless Dirac point surrounded with three 
massless Dirac points of opposite winding number, sep- 
arated by (three) saddle points at energy m*v 2 /2. A 
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Fig. 6. At low energy, the spectrum of graphene bilayer is not 
exactly quadratic, it has a fine structure made of four Dirac 
points around the each of the two points (K , —K) of the BZ. 

non-zero value of A may induce topological transitions 
of two kinds (Fig. 1-c). When A < 0, the central Dirac 
point merges with one of the three other points, leaving 
the other two unaffected. This transition is of the first 
kind mentioned above ((— 7r,7r) — > 0) and is properly 
described by Hamiltonian (5). When A > , the cen- 
tral Dirac point merges with two other points, leaving 
the fourth one unaltered [17]. Fig. 7 shows the evolu- 
tion of the relative phase <j)k and of the Berry phases 
attached to each point, in the two situations. 
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Fig. 7. Dependence of the phase <pk for the Hamiltonian (8). 

The Landau level spectrum in this model has been 
calculated recently [18]. We have given a full descrip- 
tion of this very rich spectrum, together with a semi- 
classical analysis [19]. Figure 8 shows a comparison be- 
tween the numerically obtained spectrum (full lines) 
and that (dashed lines) derived from the semiclassi- 
cal quantization rule (3) as a function of the mismatch 
parameter A. The spectrum can be divided into four 
distinct regions, according to the connectivity of the 
Fermi surface, separated by the two green curves which 
represent the energies of the saddle points connecting 
adjacent Dirac cones. The insets represent the pockets 
delimited by a constant energy in the different regions 
of the spectrum. By semiclassical quantization (3) of 
these pockets, we have been able to fit perfectly (except 
near the saddle points, where semiclassical quantiza- 
tion breaks down) the full Landaulevel spectrum. The 
different sets of levels are related to the different pock- 
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Fig. 8. Landau level spectrum for the Hamiltonian (8) as a 
function of the parameter A. In the central low energy spec- 
trum, there are w± = 4 zero energy modes corresponding to the 
four distinct massless Dirac cones. When A increases, two of 
these modes get a finite (positive and negative) energy, leaving 
only w p = 2 topologically protected modes. 



ets. Applying the quantization rule (3), it is essential 
to attribute the correct winding number wc (or Berry 
phase) to the different pockets. This is illustrated in 
Fig. 9. Finally, we find that in the central low-energy 
region, there are four zero-energy Landau levels (for 
each valley and each spin direction) corresponding to 
the four Dirac points. Among these four modes, only 
two are topologically protected, corresponding to the 
total algebraic Berry phase w p = 2 attached to these 
four points. 



Semiclassical quantization A = 2neB(n + y) 
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Fig. 9. Semiclassical quantization of the different pockets in- 
volved, with different Berry phases, eq. 3. 



5. Zero modes and semiclassical quantization 
rule 

In conclusion, we provide a heuristic argument relat- 
ing the total number of topologically protected zero- 
energy modes to the semiclassical quantization rule. 
Consider first the Hamiltonian [21,22] 



U = \ 




(9) 



describing a Dirac point with energy spectrum e = 
±A|g| p and with a winding number p (a topological 
Berry phase <J)b = =bp7r). In a magnetic field, perform- 
ing the substitution II — > \/2eB a where a is the an- 
nihilation operator, it is straightforward to find that 
the Landau level spectrum in a magnetic field is given 
by e n (B) = ±X(2eB) p/2 ^/n(n - 1) • • • (n - p + 1). In 
the limit of large n, this spectrum approximates as 
e n [B) ~ ±X[2eB(n + \ - f )] p/2 . This corresponds to 
the semiclassical quantization rule (3), wc = wt = w p 
being equal to the total winding number p. 

Consider now a Hamiltonian describing p massless 
Dirac points with Berry phase tv and p massless Dirac 
points with Berry phase — it. It can be written as (1) 

with f q = X EI -Li(n* - /?*) ]lLi(n - a*) , where on 
and fij are the complex positions in reciprocal space of 
the 7r-phase and — 7r-phase Dirac points, respectively. 
The total number of Dirac points, and thus the maxi- 
mal number of zero-energy Landau levels, is w t = p+p' '. 
In order to find the total number of topologically pro- 
tected zero-energy levels, we continuously modify the 
parameters ai,/3j —> so that fk becomes n* p IF. In 
a magnetic field, assuming for example that p > p' , 
this term keeps the form \JleB (a^a) p a p ~ p indicat- 
ing that there are w p = p—p zero energy levels. More- 
over in this case, it is straightforward to show that the 
energy levels scale as y^n(n — 1) • • • (n — p + p' — 1) so 
that the large-n expansion reproduces the semiclassical 
quantization rule with w p — p — p' . Therefore although 
the maximum number of zero-energy Landau levels is 
w t = p + p , a quantum- mechanical coupling between 
them partially lifts the degeneracy, but w p = \p — p' \ 
zero modes remain topologically protected. 
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